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1. Introduction 

Through the last decades, the theory of hyperbohc dynamical systems has become a 
cornerstone in the foundation of non-equilibrium statistical mechanics. The interest 
in hyperbolic models stems from relations between macroscopic features characteristic 
of irreversible phenomena, such as entropy production or transport coefficients, and 
dynamical properties, such as Lyapunov exponents [H [21 [3l H] . 

Billiard models, whose dynamics are conveniently described by the respective 
collision map, have proved extremely useful in this regard [5]. The simplest example is 
the two-dimensional periodic Lorentz gas with finite horizon, which exhibits diffusion 
without a drift [6]. In a statistically stationary state, this model represents a 
mechanical system at equilibrium, whose distribution, measured along the boundary of 
the scattering discs, is uniform in the position and normal velocity angles. Meanwhile the 
process of relaxation to equilibrium itself is characterized by deterministic hydrodynamic 
modes of diffusion with fractal properties [7]. One can also induce a non-equilibrium 
stationary state on a cylindrical version of this billiard by coupling it to stochastic 
reservoirs of point-particles at its boundaries [Hj; these are the so-called flux boundary 
conditions. In this difference between the chemical potentials of the reservoirs 

will transform to a linear gradient of density across the system, which is responsible for 
a steady current, given according to Pick's law of diffusion, and positivity of the entropy 
production, see [3]. 

Whereas the previous examples of billiards are Hamiltonian systems which preserve 
the phase volume, one can also consider dissipative periodic billiards driven out of 
equilibrium by the action of an external fleld. Thus the Gaussian iso-kinetic periodic 
Lorentz gas is similar to the usual periodic Lorentz gas, but for the action of a 
thermostated uniform external fleld which bends the trajectories in the direction of 
the fleld while keeping the particle's kinetic energy constant pj. This fleld therefore 
induces a non-equilibrium stationary state to which is associated a constant drift current, 
according to Ohm's law |9j. Moreover, due to the dissipation induced by the thermostat, 
phase-space volume is, on average, contracted. Therefore the sum of the Lyapunov 
exponents is strictly negative in the presence of the fleld and one can in fact identify 
this sum as minus the entropy production rate [lOj. 

Gaussian iso-kinetic dynamics are different from Hamiltonian dynamics in that the 
latter preserve the total energy, whereas the former preserve the kinetic energy. One can 
nevertheless describe the Gaussian iso-kinetic dynamics by a Hamiltonian formalism, as 
shown in [11]. A general theorem due to Wojtkowski [12] stipulates that Gaussian 
iso-kinetic trajectories are geodesic lines of the so-called torsion free connection (also 
known as the Weyl connection). The Gaussian iso-kinetic Lorentz gas can thus be 
conformally transformed into a distorted billiard table on which trajectories become 
straight lines; besides, the conformal transformation preserves the specular character of 
the collision laws [13]. Although the trajectories thus transformed do not have constant 
speed anymore, one can introduce, for any given trajectory, a time-reparametrization 
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under which the speed does remain constant. 

A volume-preserving biUiard in a distorted channel is therefore a natural 
generalization of iso-kinetic dynamics in a field-driven Lorentz channel. Ignoring the 
strict periodicity of the latter billiard, one may put aside the problem of time scales 
altogether and consider the dynamics of independent point-particles moving at constant 
speed in the new geometry. Such self-similar billiard channels were introduced in 

[15] : they consist of an infinite sequence of (nonidentical) two-dimensional cells 
that are attached together and make a (nonuniform) one- dimensional channel. The 
cells are identical in shape, but their sizes are scaled by a common factor. As a particle 
moves from one cell to a neighbouring one, its velocity remains unchanged while the 
length scales are expanded (or contracted), so that the time scales between collisions 
change accordingly. A noticeable property of these billiards concern their long-term 
statistics; even though their dynamics preserve phase-space volumes, their statistics are 
characterized by a non-equilibrium stationary state with fractal properties and a drift, 
quite similar to the Gaussian iso-kinetic Lorentz gas itself. Close to equilibrium (when 
the scaling between neighbouring cells is close to unity), one can relate the constant 
drift velocity to the diffusion coefficient of the equilibrium system, which is none but 
the usual Lorentz channel, see [T5] . 

In a recent paper [Hj, Chernov and Dolgopyat unveiled a remarkable feature of 
the drift term of self-similar billiards. Namely that the linear growth of the particle's 
average displacement with respect to time is periodic on logarithmic time scales, or 
log-periodic, with a period specified by the scaling parameter. This term may therefore 
display periodic oscillations. The purpose of this paper is to provide further insight into 
this peculiar phenomenon, which actually stems from the discrete scale invariance of self- 
similar billiards [17J. We will offer numerical evidence which demonstrates the existence 
of log-periodic drift oscillations for self-similar billiards which are sufficiently far away 
from equilibrium. Closer to equilibrium, these oscillations are, at least numerically, 
found to vanish. Thus, for all practical purposes, we may assume the drift is constant 
close to equilibrium, and grows as a linear function of the scaling parameter. Further 
away from equilibrium, this linear dependence of the average drift in the scaling exponent 
remains valid, despite the presence of drift oscillations. 

The paper is organized as follows. The model is described in section [21 which is a 
slight modification of the models previously studied in [T11[T5]. Statistical properties of 
the model are briefiy discussed in section [3l The characteristics of stationary drift, with 
numerical results, are presented in section [H Conclusions and perspectives are offered 
in section [H 

2. Description of the model 

In order to display the log-periodicity of the drift function of self-similar billiards 
discussed in section [U which was predicted in [16], we will modify the self-similar billiard 
that was defined earlier in [Ml [15], and use a unit cell rather similar to the one obtained 
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after conformal transformation of the iso- kinetic Lorentz channel [T3] . The resulting cell 
has enhanced symmetry in that all the discs are now equivalent. The main motivation 
for introducing this modification is that it allows to take much greater values of the 
scaling parameter within the constraints imposed to guarantee ergodicity. 

In order to define the self-similar channel, we start with the description of the 
reference cell, examples of which are shown in figure [H 




Figure 1. Three examples of unit cells with short trajectories. The parameter values 
(defined in the text) are respectively e = 0.5 and p = 0.459738 (left), e = 1. and 
p = 0.440862 (center), and e = 1.5 and p = 0.413317 (right). 

We define the reference cell as the interior of a region bounded on the right and left 
sides by two arc-circles, with common center and respective radii exp(±e/2)/e, and, on 
the upper and lower sides, by two oblique lines of slopes ± tan \/3e/4j|, with the exclusion 
of obstacles, which are here taken to be regular discs. The parameter e will henceforth 
be referred to as the scaling exponent, in contrast to the parameter /i which was defined 
in [m [15] (or r in [16]), and which we refer to as the scaling factor, jj, = exp(e)|§|. 

Let the origin {x, y) = (0, 0) be on the middle of the arc-circle of the left-hand side 
border. We have the central disc of radius p, located at 



I We note that these slopes may well be vertical and beyond when e > tt/VS ~ 1.814. This is not a 
restriction. 

§ The conformally transformed Gaussian iso-kinetic Lorentz channel has a similar cell geometry, but 
for the obstacles, which are deformed, slightly flattened, discs, see for details. There the scaling 
exponent is given by the amplitude of the external forcing field. 
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two discs of radii pexp(— e/2) on the inner circle, at 

(^^^^^[cos(V3e/2) - l],±^^Pi^sin(v^e/2)), 

and two discs of radii pexp(e/2) on the outer circle, at 

C-^^ cos(v^e/2) - ^^^PizfZ^, ±!^^ sin(y3e/2)) . 

Moreover, the corners where the boundaries of the cell intersect are located exactly at 
the centers of the outer discs. 

Let T>o denote the reference cell. The extended system is obtained from the reference 
cell thus defined by making a copy of it, denoted by Pi, which we scale by a factor of 
exp(e) and attach to the right-hand side of the reference cell. Likewise, another copy 
of the reference cell, this time denoted by 'D_i, is scaled by exp(— e) and attached to 
its left-hand side. By repeating this procedure, copying and scaling the right- and left- 
most cells and attaching these copies to the existing sequence of cells, we obtain the 
self-similar Lorentz channel, which consists of an infinite collection {"Pnlnez of such 
cells. Obviously the usual Lorentz channel [3] is recovered for e = 0. 

Note that when the value of p is large enough with respect to e , the central circle 
of the reference cell T>o will intersect with its inner arc-circle border. Correspondingly 
a sixth disc appears in the reference cell, whose center is located in cell T>i, to the right 
of the reference cell (as is the case in the examples shown in figure [1]). This happens 
when p > [1 — exp(— e/2)]/e. 

A point-particle which moves in this extended system has unit speed and bounces 
off elastically when it collides with the scatterers and upper and lower walls. The circular 
walls which separate the cells from one another have no incidence on the dynamics. 

As noticed in [H], the dynamics on the extended system can be boiled down to a 
single cell. Periodic boundary conditions, with rescaling of the arc-lengths and velocities, 
must then be imposed when trajectories collide on the circular sides to the right and 
left, reappearing on the opposite sides. 

3. Ergodic properties 

Although the geometry of the self-similar Lorentz channel is similar to that of the 
conformal transformation of the Gaussian iso-kinetic Lorentz gas, their dynamics are 
very different. By taking circular scatterers for all values of e, we ensure that no 
transition to a non-hyperbolic regime will occur. Nevertheless the similarities between 
the two lead to several immediate results regarding the ergodic properties of self-similar 
billiards. We refer the reader to [16] for more details. 

Two constraints have to be imposed on the parameter values, which have 
straightforward expressions in this geometry. 



Log-periodic drift oscillations in self-similar billiards 



6 



The first constraint is the finite horizon condition, by which we avoid the possibihty 
that some trajectories be baUistic. It is a straightforward generahzation of the 
corresponding condition in the periodic Lorentz gas and here becomes 

sin(y3e/4) 



P ^ Pmin — 



(1) 



The second condition is that two discs on the same radius (with center at 
exp(— e/2)/e) do not overlap so that the cell remains connected, 

1 exp(e) — 1 



P < Pmax = 



(2) 



e exp(e) + 1 

The parameter values compatible with these two conditions are shown in figure [H 




0.38 



Figure 2. Allowed values oi p v. e from equations ([T]) (p above the blue curve), and 
([2]) (p below the red curve). Ergodicity will not be guaranteed outside the domain 
bounded by these two curves. The dashed line indicates the line of parameters we use 
in the computations presented in section 2) 



4. Oscillations of the drift 

Because of the biased geometry of the system, point-particles will preferentially move 
from small cells to large cells, i. e. from left to right, thus inducing a current of mass 
along the horizontal axis of the channel. Let q{T) denote the displacement along this 
line after a time T. On average, q{T) grows linearly in T, which is to say the current 
has a steady average value in time. 

To be more precise, one can prove pTO] that the ratio q{T)/T remains of order one, 
in the sense that for any (small) p > there are constants < a < 6 < oo such that 
q{T)/T stays between a and b with probability 1 — p. However, a detailed analysis [TB] 
shows the current may actually retain a time dependence between those bounds. 

The reason is that the probability distribution of the ratio q{T)/T does not stabilize 
as T — i> oo; rather it keeps oscillating. This happens because the cell sizes grow 
exponentially; thus the size of the current cell (that is the cell where the particle is 
located at time T) is comparable to its entire previous displacement q(T). To emphasize 
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Table 1. Parameter values used in the computations whose results are displayed in 
figures [3 a and El 



e p 


e P 


e P 


e p 


0.1 0.466231 
0.2 0.465406 
0.3 0.46404 
0.4 0.462145 
0.5 0.459738 


0.6 0.456839 
0.7 0.453475 
0.8 0.44967 
0.9 0.445456 
1.0 0.440862 


1.1 0.435921 

1.2 0.430665 

1.3 0.425125 

1.4 0.419333 

1.5 0.413317 


1.6 0.407108 

1.7 0.400732 

1.8 0.394213 

1.9 0.387575 
1.99765 0.380998 



this effect, assume for a moment that the scahng factor between neighboring cells is 2 
(z. e. the scaling exponent e = In 2). Then the current cell would be just as big as all 
the previous cells combined. In other words, the geometric shape of the current cell 
scales up at the same rate as the particle's entire path. 

This observation should convince the reader that the geometric structure of the 
current cell might affect the macroscopic evolution of q{T)/T. For example, there may 
be some parts of the cell which, on average, the particle traverses faster than others; 
in that case, as it passes through 'fast terrain', the ratio q{T)/T would grow, and as it 
drags itself through 'slow motion areas', the ratio q(T)/T would decrease. This happens 
periodically, as the particle moves from one cell to the next (even bigger) cell. 

Note that the period of these oscillations is not fixed; it keeps growing exponentially 
with T. More precisely, the length of the period corresponds to the time it takes the 
particle to traverse the current cell completely, thus the ratio logg(T)/e should have 
period one. 

These arguments were made precise in [16], where the following theorem was 
proved : 

Theorem 1 Assume e is chosen such that equations (QP and ^ are satisfied. The 
following holds. Let {T„}„gM be a sequence of increasing times, Tn oo, such that 
the fractional part o/logT^/e tends to a constant < 5 < 1. Then the distribution of 
q{Tn)/Tn converges to a limit. 

The theorem also implies that the limit of q{Tn)/Tn may depend on the limiting 
fractional part of logT„/e, i. e. of 6. We may therefore expect that {q(T))/T will display 
oscillations as a function of logT/e, with unit period. On the other hand, there may be 
no oscillations. Whether we can observe this oscillatory regime or not depends on the 
parameter values, as our numerical analysis reveals. 

For the purpose of numerically demonstrating the drift oscillations, we let the 
scaling exponent e to vary between and its maximal allowed value near e = 2 (see 
figure [2]) and take the corresponding radius to be p = l/2(pmax + Pmin), in between the 
two bounds specified by equations ([H) and (I2l). The precise values of the parameters we 
used in our computations are shown in Table [TJ 
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Figure [3] shows the results of numerical computations of the drift function {q{T)) /T 
using 10^ trajectories with random initial conditions located on the central circle (or on 
the discs on the horizontal line when there is an overlap) of cell Vq, and unit velocity at 
random angle. For each of these trajectories, we computed the horizontal position q{T) 
at times T which we took to be logarithmically spaced on the interval < logT/e < 100, 
so as to have 25 points on each unit interval of the scale logT/e. 




Figure 3. Average drift terms {q{T))/T vs. log(T)/e. Here shown are the parameter 
values of Table [1] corresponding to e = 0.8, 0.9, 1.0, 1.2, 1.4, 1.6, 1.8, 1.9 et 1.99765, 
in order of increasing e, left to right and top to bottom. The dashed lines show the 
computed averages, as plotted in figure [5] below. The color code is a hue, regularly 
varying from cyan to red as units of logT/e increase from 50 to 100. Note that the 
span of the a;-axis was reduced for the larger e values in order to better display the 
periodicity of the oscillations. 

For values of the scaling exponent e > 1.0, the regularity of the observed periodic 
oscillations is spectacular, as confirmed by the collapsed curves of figure HI On the 
contrary, for values of e < 0.8 (only e = 0.8 is shown on figures [3] and HI), our 
measurements do not reveal periodic oscillations predicted by theorem [1], which, in 
principle, applies to all values of e. The lack of noticeable periodic oscillations for small 
e is, most likely, the effect of insufficient time period T in our simulation (recall that 
theorem [1] only holds asymptotically, as T — >^ oo). It is (theoretically) possible that, after 
many more (say, 10^ or 10^) initial periods we would see proper periodic oscillations. 

The summary of our computational results is shown in figure O The left-hand panel 
shows the average drift as a function of e, computed by taking the average of the time 
series of {q(T))/T for every parameter value of Tabled! The right-hand panel shows 
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Figure 4. Same as figure [3] with the horizontal axis taken modulo 1. For every figure, 
we plot the last 50 oscillations. The color code is identical to figure [3l When e is large 
enough, the curves nicely collapse on the same periodic function. 



the computed average amplitude of the oscillations, which we denote by A, with the 
corresponding error bars. The amplitudes of the oscillations is here measured by taking 
the average of one half of the sum of the maximum and minimum values of {q{T)) /T 
over each unit interval of logT/e, which corresponds to a period of oscillation. Values 
of e < 0.7 were discarded, because oscillations could not be observed (yet). 



{q(T)) 
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Figure 5. Average drifts (left) and oscillation amplitudes (right) measured for the 
parameter values of table [TJ 

The latter graph clearly demonstrates a power law behavior of the amplitude in 
the scaling factor, i. e. A ~ exp(ae), with the coefficient a ~ 1.5, for parameter values 
1.5 < e < 2.0. The coefficient seems to be somewhat larger for smaller values of e. The 
reason for this difference we believe is due to a non-trivial dependence in the parameter 
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p, which here varies non-hnearly with e. It might also be due to our hmited statistics 
(as also reflected by the large error bars). 

5. Concluding remarks 

To conclude, self-similar billiard channels are instances of non-equilibrium chaotic 
billiards with volume-preserving dynamics, whereby a geometric constraint induces a 
current of mass, going from the smaller to the larger scales. The paper clearly establishes 
the following three novel results : 

(i) The average drift is a linear function in the scaling exponent e throughout its range 
of allowed values. 

(ii) When the scaling exponent is large enough, log-periodic drift oscillations do occur, 
i. e. the drift is not a constant. 

(iii) The amplitude of the oscillations follows a power law in the scaling factor, more 
clearly so for the larger range of parameter values. 

It is to be noted that property (ii) is an expected consequence of the discrete scale 
invariance of the system, which typically manifests itself in the presence of power laws 
with complex exponents. The signature of such power laws is the log-periodic corrections 
to their scalings, in this case the corrections to q{T) ~ T. Similar phenomena are 
observed in many different physical situations, see [17] and references therein. That 
the drift oscillations are barely noticeable when the scaling exponent decreases below 
e = 1.0 is perhaps not surprising as log-periodicity is typically a very small effect. Going 
to the other end of the parameter range, it is actually remarkable that the regularity of 
drift oscillations is so well pronounced for larger values of e. In this regard, self-similar 
billiards offer new ground to further studies of this interesting phenomenon, here in the 
framework of hyperbolic dynamical systems, and thus improve our understanding of a 
ubiquitous property which has already found many known flelds of applications. 
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